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Abstract 

We study the LHC signal of a light gluino whose cascade decay is dominated by channels 
involving top, and, sometimes, bottom quarks. This is a generic signature for a number of 
supersymmetry breaking scenarios considered recently, where the squarks are heavier than 
gauginos. Third generation final states generically dominate since third generation squarks are 
typically somewhat lighter in these models. At the LHC we demonstrate that early discovery 
is possible due to the existence of multi-lepton multi-bottom final states which have fairly low 
Standard Model background. We find that the best discovery channel is 'same sign dilepton'. 
The relative decay branching ratios into tt, tb and bb states carry important information about 
the underlying model. Although reconstruction will yield evidence for the existence of top 
quarks in the event, we demonstrate that identifying multiple top quarks suffers from low 
efficiency and large combinatorial background, due to the large number of final state particles. 
We propose a fitting method which takes advantage of excesses in a large number of channels. 
We demonstrate such a method will allow us to extract information about decay branching 
ratios with moderate integrated luminosities. In addition, the method also gives an upper 
bound on the gluino production cross section and an estimate of the gluino mass. 
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1 Introduction 



Many considerations of new physics at the TeV scale point toward top rich final states at the 
LHC. Such scenarios include top compositeness [JJ, and models in which top partners ensure the 
naturalness of electroweak symmetry breaking [5], [31 HI [H] • Compared with the Standard Model 
QCD production of tt, top quarks in the final states of new physics production typically either have 
very different kinematics [51 [7J [5] , or different event topology [HI El HH QJ2 H3] j which may make it 
crucial to develop new techniques to identify them. 

Naturalness of the electroweak symmetry breaking typically requires the existence of a light 
top partner. Experimental observations, ranging from the existence of Cold Dark Matter in the 
universe to the constraints of electroweak precision measurements, strongly motivates the existence 
of a neutral stable particle as particle of new physics signals at the LHC. The above considerations 
lead to the recent studies of new physics signals with ti + ]£t final state [9] • 

In this article, we consider scenarios which also include a gluon partner, such as gluino in low 
energy supersymmetry [31 [T3], the KK-gluon in universal extra dimension models |3J, or other 
octet states [T5J [TH1 U7J . Due to the nature of proton colliders, production of such color octet gluon 
partners and their decay typically becomes the main channel of new physics signals. Decay products 
of the gluon partner typically include an even number of quarks. Combining this with the scenario 
with light top partners, we conclude that a typical signature of production of gluon partners will 
be multiple top quarks in the final states. 

In the rest of this paper, we will study this signature using a particular example, low energy 
supersymmetry with light gluinos. We will focus on the scenario in which the squarks are heavier 
than the gluino and the third generation squarks are lighter than those of the first two generations. 
In this case, the gluino will dominantly decay into top (bottom) quarks. Although not absolutely 
unavoidable, this is a generic possibility from the point of view of many models being studied 
recently [HI [TBI E01 EJJ ■ Heavier squark masses are often preferred due to constraints from flavor 
changing neutral currents [22j E31 [24] . RGE running of scalar masses from the high scale down to 
the electroweak scale will tend to push the third generation squark masses significantly lower than 
those of the other generations. Large third generation trilinear couplings will also help further lower 
one of the stop masses. For earlier studies on gluino decay into third generation quarks, see Refs. 

1231251 EE 123 ESI EDI EH- 

In the minimal case, this scenario has only a single, light top partner, tn. On the other hand, 
we would like to include in our study the possibility that all the squarks of the full third generation 
could be light relative to the other two generations. Therefore, we are led to consider decay channels 
g — > tiN , g — > tbC~ , and g — > bbN. As will be clear from our discussion later, many of the leading 
order features of the signature of different combinations of the above decay channels (from gluino 
pair production and decay) can be quite similar. Therefore, it is one of the primary purposes of 
this article to study techniques to distinguish them. We remark that such measurement is crucial 
for understanding both the spectrum of the third generation sfermions and the electroweak-inos. 
For example, if we measure a non-zero branching ratio for the decay channel g — > tb\~ , then there 
must be a light electroweak-ino carrying charge, suggesting that the lightest supersymmetric particle 
(LSP) is mostly Wino. Moreover, the relative branching ratio BR(g — > tt + \)/BR{g — ► bb + x) 
carries important information about the squark masses. For example, if this ratio is close to 
1, it strongly suggests that left-handed squark masses are lighter than the right-handed ones, 

m Q 3 < "V ^ 

In section El we will focus on discovering new physics in this class of final states. Decay of 
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multiple top quarks could lead to b-rich and lepton rich final states. Therefore, we expect great 
potential for early discovery. For example, we show that significant excesses can be observed in 
many channels even with just 500 pb _1 of data. The obvious channel with the best early discovery 
potential is same-sign dilepton plus additional b-tags. 

In section 21 we study the problem of direct reconstruction of top quarks. Large combinatorics 
and high probability of object merging are expected due to the large multiplicity of final state 
particles. Therefore, while we may gather evidence from this study that there are some top quarks 
in the decay chain of the gluino, similar to the approaches taken in Ref. [T5] , the direct 

reconstruction efficiency for the top quark is low. We find that while the efficiency for reconstructing 
a single top quark candidate approaches 49%, the efficiency to detect three and four candidates 
drops dramatically to 1.5% and 0.02%, respectively. As a result, it is less likely we can measure top 
multiplicity by direct reconstruction. 

We will demonstrate in our study a fitting procedure which could allow us to measure the 
branching ratios of different gluino decay channels. We begin by simulating a number of samples of 
gluino pair production and decay, each with different final states, such as tttt, tttb, ttbb, and so on. 
Then we will fit the relative weights of different samples to match a set of experimental signatures. 
Of course, without precise knowledge of the underlying spectrum, choice of the templates will 
introduce errors in the estimate of the branching ratios. We studied such effects by using several 
templates with different hypotheses for the relevant masses. We conclude that such a method will 
allow us to establish important features of gluino branching ratios. 

We carry out our study on several benchmark models with relatively low gluino masses. A 
detailed scan of the parameter space involving the gluino mass and different branching ratios is 
beyond the scope of this paper. The corresponding results for heavier gluino masses (but with 
similar decay branching ratio and mass difference between gluino and the LSP) could be roughly 
obtained by scaling from the present result using relative production cross sections. The mass gap 
between the gluino and the neutralino or chargino in the next step of the decay chain could also 
have important effects as it will affect the detection efficiency of various decay products. In general, 
a larger mass gap will enhance the discovery potential. At the same time, we expect this effect is 
milder in comparison with the dependence on the gluino mass. 

We emphasize that our goal in this study is to demonstrate a method which allows us to extract 
information of the SUSY spectrum, such as the identity of the LSP with relatively low integrated 
luminosity. This is possible mainly because, unlike the precision measurements of the masses and 
couplings, our method mainly relies on inclusive counts and general kincmatical features. Moreover, 
since we do not demand direct reconstruction, we are able to take advantage of many channels with 
multiple leptons. Due to lower background and theoretical uncertainty in comparison with the pure 
hadronic channel, we expect to have significant excesses in many of these channels. After discovery, 
we expect our method will yield a first set of clues about the underlying model during the early 
stages of LHC operation. 

2 Benchmark Models 

We consider four benchmark models in our study. The model parameters and relevant decay 
branching ratios are shown in Table [T] For simplicity we will only consider gluino decay chains in 
which the only observable decay products are tops and bottoms, and they only come from one single 
step in the gluino decay chain. This is of course the case if the only state lighter than the gluino 
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Model parameters (TeV) 


Branching ratios 




m~ a 


m «3l,2 


m h 






m b 2 


m N,C 


(g - tt") 


(.9 -> W.) 


(g - tt) 


A 


0.65 


8 


1.3 


8 


2.5 


8.1 


0.1 


0.92 


0.07 





B 


0.65 


4 


0.8 


0.93 


0.87 


4 


0.1 


0.71 


0.27 





C 


0.65 


4 


0.64 


0.9 


0.72 


4 


0.1 


0.52 


0.47 





D 


0.65 


4 


0.63 


0.9 


0.72 


4 


0.1 


0.09 


0.22 


0.69 



Table 1: Model parameters and relevant branching ratios for the benchmark models considered in 
this paper. The mass parameters are in TeV. The models A, B, and C have bino LSP. In Model 
D, the lightest neutralino and lightest chargino are both winos. We have adopted the short hand 
notation where we omit the explicit mention of the identity of the electroweak ino in the decay, as 
it can always be inferred from the observable particle content. 



is the neutral LSP, i.e., a bino-like LSP (Models A, B and C). We also consider the possibility of a 
wino LSP, in which the neutral LSP is almost degenerate with the lightest chargino (Model D). 

Model A is the simplest example of multi-top physics. It is designed to have only a single light 
stop, and therefore always produces four tops in the final state. Model B and C are designed to 
include the decay channel g — ► bbN, with each model exhibiting a somewhat different branching 
ratio. In Model D, the wino LSP is approximately degenerate with the lightest chargino, which 
is also wino-like. It is designed to further include a chargino in the decay chain. We will study 
distinguishing it from the other channels. Since the charged wino is approximately degenerate with 
the wino LSP, it appears only as missing energy. 

In Table [1] and the rest of this paper, we will adopt a short-hand notation for gluino decay by 
only including top and bottom quarks and not giving explicitly the electroweak-inos, as it should 
be evident from the context. 

3 Signal Isolation and Backgrounds 

We begin by discussing the prospects for signal isolation above Standard Model background at 
the LHC. The relatively large 6-jet and lepton multiplicity associated with four-top production 
provide for potentially striking signatures that are easily distinguishable above the expected SM 
background. We find that by requesting > 3 6-tagged jets and at least one lepton, it is possible 
to achieve signal significance S/VB > 3 for only 500 pb _1 of integrated luminosity. We will 
demonstrate the discovery potential in three of the four benchmark models. 

One of the important backgrounds from the Standard Model for final states with many 6-tagged 
jets, several isolated leptons and very high missing Et, is top pair production, tt. The expected 
cross-section at the LHC for this background is u — 833 pb (NLO+NLL result 32]). The tt 
background event samples were produced using Pythia 6.4 [33] . 

We have also included in our analysis a set of SM backgrounds involving associated production 
of W/Z bosons with third generation quarks. These contribute significantly to signals with high 
lepton multiplicity, or same sign dileptons in the final state. As we will see, the latter case is a 
particularly important discovery channel early on. The parton-lcvcl SM background event samples 
were produced with Madgraph v. 4. 2. 3 [34], with the exception of the tt background which was 
produced using Pythia 6.4. The tt cross section was taken from [32], while the cross sections 
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Process 


a[pb] 


Process 


a[pb] 


Process 




tt+ 1,2,3 jets 


833 


tbZ + 1,2 jets 


0.67 


ZW + b + 1,2 jets 


0.48 


ttZ+ 1,2 jets 


0.28 


+ 1,2 jets 


0.58 


1,2 jets 


0.50 


UW~ + 1,2 jets 


1.5 


t6VF+ + 1,2 jets 


0.18 


+ 1,2 jets 


0.85 


tt¥+ + 1,2 jets 


3.4 


tbW~ + 1,2 jets 


0.09 


+ 1,2 jets 


0.28 



Table 2: Standard Model backgrounds and relevant cross sections used. 



Standard Model Background 





0b 


lb 


2b 


3b 


>4b 


0L 


1717.46 


3069.29 


2091.27 


320.54 


36.51 


1L 


783.34 


1489.85 


998.8 


118.42 


8.49 


OS 


41.89 


61.82 


34.06 


4.46 


0.01 


ss 


0.41 


0.97 


0.44 


0.04 




3L 


0.1 


0.54 


0.24 


0.06 





>4L 



Table 3: Number of Standard Model events with n (n = 0..4) b-tagged jets and m (m = 
0, 1, OS, SS, 3, 4) leptons for the combined SM background considered. The following cuts were 
applied: MET > 100 GeV, at least 4 jets with p T > 50 GeV, all jet and lepton p T > 20 GeV. The 
results are normalized to 500 pb^ 1 . 

for the other backgrounds are calculated from Madgraph. The subsequent parton shower and 
hadronization were simulated by Pythia 6.4. We have used the CKKW matching scheme [35] 
implemented in Madgraph. The events are then passed on to PGS-4 [35] with parameters chosen 
to mimic a generic ATLAS/CMS type detector. All background sources considered, and their 
respective cross sections are given in Table [2] 

The signal event samples, for gluino pair production and decay, were produced using Pythia 
6.4 and have been passed through the same PGS-4 detector simulation. Appropriate k-factors [37] 
were applied to the LO signal cross-section calculated by Pythia to obtain the NLO cross-section. 

Basic muon isolation was applied to all samples: If the summed Pt in a AR — 0.4 cone around 
the muon is greater than 5 GeV, or the ratio of the Et in a 3 x 3 cell region of the calorimeter to 
the muon Pt is greater than 0.1125, the muon is merged with the nearest jet in AR. 

We have also imposed on both the signal and the background the following selection cuts 

1. E T > 100 GeV 

2. pt > 20 GeV and pseudorapidity |?7| < 2.5 for all objects 

3. At least 4 jets with p T > 50 GeV 

Table [3] shows the expected number of events from the SM background. We have classified them 
according to the number of 6-tagged jets and isolated leptons in the event. Same sign (SS) and 
opposite sign (OS) di-leptons are separated as they have very different origins and sizes. We will use 
the possible excess in these channels to assess the discovery potential. The results are normalized 
to 500 pb^ 1 of integrated luminosity. Crossed out entries indicate no background events passing 
the signature and selection cuts. 



5 



Model A Model C Model D 





2b 


3b 


>4b 




2b 


3b 


>4b 




2b 


3b 


>4b 


1L 




166.2 


70.9 


1L 




106.5 


44.2 


1L 




98. 


37.8 


OS 


27.6 


19.3 


7.3 


OS 


13.3 


10. 


3.9 


OS 


5.6 


3.8 


1.5 


ss 


12.7 


9.4 


3.1 


SS 


4.1 


2.8 


1. 


SS 


2.9 


2.2 


0.8 


3L 


3.1 


2.2 


0.8 


3L 


1.1 


0.6 


0.2 


3L 


0.2 


0.1 


0.1 



Table 4: Number of signal events passing the selection cuts and containing n b-tagged jets and m 
leptons. The selection cuts applied were: MET > 100 GeV, at least 4 jets with pt > 50 GeV, all 
jet and lepton pt > 20 GeV. The results are normalized to 500 





Model A 






Model C 






Model D 






2b 


3b 


>4b 




2b 


3b 


>4b 




2b 


3b 


>4b 


1L 




15.3 


24.3 


1L 




9.79 


15.2 


1L 




9.00 


13.0 


OS 


4.73 


9.12 


87.0 


OS 


2.28 


4.73 


47.0 


OS 


0.957 


1.79 


18.3 


SS 


19.2 


49.4 




SS 


6.10 


14.5 




SS 


4.31 


11.3 




3L 


6.44 


9.26 




3L 


2.35 


2.63 




3L 


0.418 


0.318 





Table 5: Signal significance S/yB, computed for the results in tables H] and [3] The crossed out 
entries indicate no background events passing the signature and selection cuts. 

Table [4] shows the expected number of signal events with n 6-tagged jets and m isolated leptons 
(leptons = e^,^). Model A, which is predominantly a four top signal, has significantly more 
muti-lepton and 6-jet events passing selection cuts than Model C and Model D, which have fewer 
four top events. Model C is a stronger signal than Model D, which has very few four top events. 

In Table [SJ we show the signal significance achievable with 500 pb^ 1 integrated luminosity. By 
requesting > 3 fr-tagged jets it is possible to observe signal significance S/vB >3 for events with 
multiple leptons, an excess consistent with multi-top production. In the single lepton channels, 
a more detailed study of the background would be required to carefully calculate the expected 
significance, since there are likely to be significant background contributions from QCD processes 
with a faked lepton. The same-sign dilepton channel is probably the best channel for discovery, a 
finding that is consistent with results in [38] . It can also be observed that already with 100 pb _1 
integrated luminosity a 4 top signal (Model A) may be established in the same sign dilepton, 3 6-jet 
channel. 

During the early period of LHC data taking, missing energy may not be well understood since it 
requires a 'global' understanding of the ATLAS/CMS detectors. Therefore missing energy should 
not necessarily be taken as a reliable tool to discover new physics at low luminosities. If we do not 
include the missing energy cut in our analysis, then QCD backgrounds, particularly bb production, 
becomes a significant background to the multi-top signal. Requiring four hard jets, as we have 
done here, does reduce the QCD backgrounds since hard jets are less likely to produce isolated 
leptons [39]. However, even though the 2-lepton background from QCD might still be significant, 
the 3-lepton QCD background will unlikely be more significant than the 4-top signal 39, 40J. Thus, 
it seems reasonable that, without using missing energy, discovery could still be possible without 
missing energy at integrated luminosities greater than or equal to 500 pb _1 . 
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4 Direct Reconstruction 



Once evidence is obtained for an excess beyond the Standard Model in events with multiple 6-jets 
and leptons, it is natural to assume that the signal involves production of multiple top quarks. 
In order to provide concrete evidence for this, W bosons and top quarks in the signal should be 
reconstructed. 

In the model under consideration, where each signal event has four top quarks, the main sources 
of difficulty in reconstruction are low statistics, large combinatorial background, and poor object 
reconstruction due to the extremely complex event topology. Every tttt event has four W bosons, 
each of which gives two jets if it decays hadronically. Furthermore, every top decay itself gives a 
6-jet. Therefore the expected number of hard jets arising in a decay with k W bosons decaying 
hadronically is 4 + 2k. On top of this there are also jets arising from initial/final state radiation 
(ISR/FSR). For comparison, ft all hadronic decays are expected to have 6 hard jets on average 
before ISR/FSR effects are included. 

We expect that due to the large number of particles in the final state, the chance of recon- 
structible objects (such as jets) overlapping is large, so that the detector will frequently be unable 
to accurately reconstruct isolated objects. If two partons are very close in AR (roughly less than 
twice the diameter of the cone used in the jet algorithm), the jets coming from the partons are 
likely to be recombined into a single jet. Similarly, if a lepton is very close to a jet, it will likely not 
pass the isolation requirement. We utilize the standard cone jet reconstruction algorithm imple- 
mented in the PGS simulator, however we reduce the AR cone size to 0.4 in line with the expected 
performance of the ATLAS detector. 

Figure Q] shows the distributions of the lowest, 2nd lowest, 3rd lowest, and 4th lowest AR values 
between reconstructible partons in four top events from benchmark model A. For comparison, the 
same information is also shown for Standard Model tt events. Here we define a reconstructible 
parton as any first /second generation quark or 6-quark. It may be seen in figure [1] (a) that in four 
top events there is a large likelihood of three to four pairs of overlapping objects, rendering reliable 
final state reconstruction difficult. 

Direct reconstruction of top quarks in somewhat different decay chains has been studied in 
Ref. [26l [27]. We expect a similar study in our case will also yield at least some evidence that 
there are indeed top quarks in the event. Rather then presenting a detailed analysis here, we focus 
here on a somewhat different question. Clearly, in order to completely measure the branching 
ratios into different final states involving different numbers of top and bottom, it is not enough to 
just reconstruct a certain top quark. We need to be able to reconstruct all of the top quarks in 
the event with reasonable efficiency. However, our study already shows that there is a significant 
overlap between different object in this type of signal events. Including additional combinatorics, 
we expect a very low efficiency for reconstructing multiple tops. 

To gain some estimate of such efficiencies, we study how many tops we can possibly reconstruct 
in a event. To proceed, we define a 'top candidate' to be the combination of two light-jets (ie non- 
6-tagged) and one 6-tagged jet, where the non-6-tagged dijet invariant mass satisfies 65 < rrijj < 95 
GeV, the invariant mass of the b-jet with any lepton must satisfy mu > 155 GeV, while the invariant 
mass with either of the two non-6-tagged jets must satisfy rribj < 160 GeV. Finally, the invariant 
mass of the final three-jet combination must satisfy 125 < rrijjb < 225 GeV. 

Figure [2] shows two distributions for the number of top candidates observed in our benchmark 
model A. The first figure, [2] (a), includes all possible three-quark combinations that satisfy the 
requirements above. The inherent combinatorical background due to the intense hadronic activity 
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(a) (b) 

Figure 1: Distributions of the lowest four values of separations, AR, between reconstructible partons 
(which here we take as any first/second generation quark or b-quark) obtained using the Pythia 
parton-level (truth) information for benchmark model A. For each event in the simulation, AR was 
computed for all pairs of reconstructible objects, ranked, and the lowest four values binned into 
histograms. The resulting distributions showing the lowest, 2nd, 3rd, and 4th lowest AR values 
encountered are given for (a) four-top events of model A, and also (b) hadronic tt decays. All 
distributions are normalized to unity. The partons present in four-top events are significantly closer 
to one another relative to those in SM tt events, increasing the liklihood of overlap, as well as a 
lower reconstruction efficiency. 

in four-top events is clearly visible. We see that the same &-quark can be combined with other 
partons to form several "top candidates" . There can be no more than 4 top quarks present in the 
event, and the vast majority of 'candiate' combinations are incorretly chosen. Figure [2] (b) shows 
the same information except that here we isolate only the distinct jet combinations of each event. 
Degeneracies that arise are removed by keeping track of the mass difference rrijjb — rn t for each 
jet triplet, and choosing the set of triplets with the lowest average difference. From the figure, it 
is clear that this approach gives a significantly more reasonable result. However, notice that the 
number of reconstructed candidates drops dramatically as the top multiplicity increases, rendering 
statistical analysis essentially impossible without a large integrated luminosity. 

The study we perform here is not an actual reconstruction of the top quark. A true reconstruc- 
tion involves positive identification of the top quark through statistical determination of the top 
invariant mass. Instead, this study is an estimate of how many objects obtained from recombin- 
ing final states can be consistent with a top quark. We expect this study, though not completely 
precise, does capture the main effect of combinatorics and object merging. We observe that the 
efficiency for detecting one top quark as defined is approximately 48.5%, for two quarks ~ 16.7%, 
three quarks ~ 1.5%, and for four quarks ~ 0.02%. 

Further perfections of reconstruction techniques are certainly going to improve these results and 
should be pursued. However, it is likely that to extract statistically significant information based 
on a full reconstruction will still require large integrated luminosity. 
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(b) 



Figure 2: As a demonstration of the combinatoric background, in figure (a), we show the observed 
number of possible top 'candidates' obtained for our benchmark model A. (a) We consider com- 
binations with one 6-jet and two light jets, where the invariant mass falls within the top mass 
window, and where the combination satisfies a minimal set of selection criteria (see text). Due to 
the combinatorics, the same 6-quark can be combined with other partons to form multiple "top 
candidates" . In (b) we show the resulting number of top candidates after attempts are made to 
remove combinatorics by isolating distinct 3-jet combinations. For this basic study we only require 
that events have at least 4 6-tagged jets. 

5 Understanding multi-top final states 

As we have demonstrated in the previous sections, new physics signals containing multiple tops 
can be discovered at the LHC. By direct reconstruction, we expect to gather evidence that there 
are indeed top quarks in the signal. However, there are many possible event topologies which 
can contribute to our signal. In general, gluino can decay into third generation quarks through 
g — > tt + Ni, g — > bb + Ni, and g — ► tb + C~ . Therefore, final states coming from gluino pair 
production can involve from zero to four top quarks, with relative amounts determined by gluino 
branching ratios. As we have argued in the introduction, measuring such branching ratios plays a 
central role in understanding the properties of superpartners involved. For example, the relative 
ratio of g — > tt and g — > bb could give us important information about the spectrum of the third 
generation squarks. At the same time, significant decay branching ratio of g — > tb strongly suggests 
that either Higgsino or Wino (or both) is lighter than the gluino. 

As demonstrated in the previous section, measuring the branching ratios by directly reconstruct- 
ing top quarks suffers from low efficiencies. In this section, we will instead tackle this question from 
a different approach. We assume that the new physics signal has already been discovered in a set 
of channels, particularly those with multiple leptons and multiple bottom quarks (Table [4}. In 
addition, there are a set of statistically significant experimental observables defined on the set of 
signal events. This general approach can in principle be applied in practice with an arbitrary set of 
experimental observables based on availability. Here, we will only use a limited set of experimental 
variables 



• 2 6-jets and either OS di-leptons, SS di-leptons, 3 leptons, 4 or more leptons 



• 3 6-jets and either 1 lepton, OS di-lcptons, SS di-leptons, 3 leptons, 4 or more leptons 



• 4 or more 6-jets and either 1 lepton, OS di-leptons, SS di-leptons, 3 leptons, 4 or more leptons 

We will consider a general set of decay channels of the gluino which can in principle contribute 
to the new physics signal. We will consider as free parameters the relative branching ratios of those 
channels whose values are determined by a fit to the above set of experimental observables. This 
approach can be viewed as a natural application of the method proposed in Ref. [43j . The set of 
possible decay channels are chosen as follows 

99 -> ttxtt_ gg -> ti_xtb_+ c.c gg -> tfybb 

gg — > tb\bb + c.c. gg — » bbxbb gg — ► tbyib + c.c. gg — > tbxtb + c.c. 

where we continue to used the short hand notation of not explicitly displaying either the lightest 
neutralino, or the lightest chargino. We have assumed that the lightest chargino and neutralino are 
nearly degenerate so that transitions between them will not yield observable decay products. +c.c. 
indicates that we included the charge conjugated event. 

The number of events in the possible gg decay channels are given by 

n t iti = <?gg£Br(g -> tt)Br(g -> ft) 



CBr(g -> tt)Br(g -> tb) 



tub — u gg 



n qqqq = a gg CBr(g -> qq)Br(g ->■ qq). 

These number can then be used to estimate the number of observed events with a particular 
signature, N Ji, which can receive contributions from several channels listed above with a particular 
fraction, e Channel' depending on event topology and experimental efficiencies. We have 

N 0b0l _ _ ObOl , _ ObOl , , OfcOJ 

Jv obs — 'HtttHttt T 'Httb^titb T • • • T ll qqqqt qqqq 

N lbOl _ _ lbOl , „ _ _,160j _j_ , AbOl 

Jv obs — 'HtttHttt ~r 'HttbHitb ' ' ' " , qqqq t qqqq 



N 4b4l _ _ 4b4l , _ 464; , , 4&4Z 

Jv obs — 'HtttHttt ' ll tttb t tttb T • • • -T li qqqq t qqqq . 

In the rest of this section, we will first obtain estimates of all signal efficiencies Charmer Then, we 
will perform a x 2 -fit to determine a set of best fit values of y^a^C -Br c hannci- Note that in order 
to obtain the branching ratio from these counting signatures, we will have to know the product 
(JggC. Such information could be available independently from other measurements, such as the 
gluino mass. We will show later in this section that our method indeed gives us an estimate of this 
absolute rate. At this moment, we note that we can already derive a lot of information about the 
underlying model if we measure the ratio of branching ratios, for which the dependence on (TggC 
drops out. 

The key step in this method of fitting for branching ratios is, of course, to obtain an estimate 
the efficiencies, channel- Many factors enter in such an estimate. The existence of many objects, 
in particular jets, in the event means that simple extrapolation from single object efficiency is no 
longer reliable. In particular, there is a significant chance now for the leptons from top decay to 
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Figure 3: Results for fits of ratio Br(<7 — * tt)/Br(g — > bb) for benchmark Model-B. The left panel 
shows the result with the correct mass hypothesis. In the right panel, the efficiencies are calculated 
for five different mass templates show in the legend. The solid horizontal line gives the actual values 
of the branching ratios. Errors are la and include errors for subtracting off the tt background. 



be close to a jet and therefore fail the isolation cut. 6-tagging can also be affected. Numerical 
simulation using an appropriate model, a template, of the underlying physics involved is therefore 
unavoidable. 

In practice, certain assumptions about the underlying model must be made in choosing a tem- 
plate. We have already chosen the set of channels to include, based on information gained from the 
type of study in the previous section. In addition, we have to choose the spectrum, the gluino mass 
and LSP mass, to be used in the template. To begin with, we will first choose a template model 
using the actual gluino mass and LSP mass as the underlying model. We will fit the underlying 
branching ratios by using efficiencies obtained by simulating this template and demonstrate such a 
fit does give us accurate information of the underlying model. We will then simulate a set of different 
templates with different mass hypotheses. We will show that for the variation of mass hypotheses 
we have studied, the difference induced for the fit does not significantly affect our conclusion about 
the underlying models. Of course, such a result still leaves the question of whether the range of 
variation of mass parameters we have used is too optimistic. Later in this section, we will show 
that indeed it is reasonable to expect we can get estimates of mass scales from other experimental 
obscrvables so that we can choose our mass hypothesis with an error within this range. For the 
continuity of our presentation, we will present our result for the fits first. Then we will come back 
to address the question of choosing the mass hypothesis in detail. 

Our fit to Br(<7 — > tt)/Br(g — > bb) for benchmark model-B using the correct mass hypothesis is 
shown in the left panel of Fig. [3l For this and all the other fits presented in this section, we have 
required that for a channel to be included in our fits, there must be at least 1 signal event, and the 
significance over the tt background must be greater than 3. 

We pause here to briefly describe how error bars are calculated in these and the other fits 
presented in this section. Using large statistics for the template model, the statistical errors for 
calculating the efficiencies and for determining the number of events for a given signature are 
assumed to be much smaller that the expected Gaussian errors from the LHC data. The only 
errors we include here are statistical errors from the minimization procedure. Therefore, the lcr 
error for a given branching ratio is given by the change in the branching ratio required to shift the 
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Figure 4: Results for fits of ratio Br(g — * tt)/Bi(g — > bb) for benchmark Model-C. The left panel 
shows the result with the correct mass hypothesis. In the right panel, the efficiencies are calculated 
for five different mass templates shown in the legend. The solid horizontal line gives the actual values 
of the branching ratios. Errors are la and include errors for subtracting off the it background. 
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From the left Fig. [3l we see that using the correct mass hypothesis, we will be able to measure 
the ratio Br(g — > tt)/~Br(g — > bb) with good accuracy for an integrated luminosity of ~ 5 — 10 fb _1 . 
In particular, we will be able to verify that in Model B, the gluino decay is dominated by g — ► tt 
with a smaller but non-vanishing branching ratio for g — * bb. 

Next, we want to assess the effect of changing our assumptions of underlying spectrum. We will 
assume that although the mass spectrum cannot be precisely measured, some crude estimates can 
still be made based on kinematical variables such as M c g and rate. We will provide justifications 
for this assumption later in this section. Therefore, we will consider cases where the gluino mass 
only deviates from the underlying benchmark model by about 100 GeV. In particular, we use four 
additional sets of alternative templates and carry out the same fit. The result in Model B is shown 
in the right panel of Fig. [3] We see that using different mass hypotheses does make a visible 
difference. However, we observe that these differences arc not big enough to dramatically affect the 
information we will extract from our measurement of Br (5 — * tt)/Br(g — ► bb). In the case of 3-body 
gluino decay under consideration, this ratio is proportional to (mi/wj) . Therefore, we see that 
using different mass hypotheses within this range will at most result in a factor two error in the 
measurement of the ratio of the branching ratios, will induce at most ~ 20% shift in the inferred 
ratio mg/my. 

Notice that incorrect assumptions about the underlying spectrum do not lead to significant 
effects in the fitting when considering at least 1 fb _1 integrated luminosity of data. In all cases, we 
can still extract a good estimate of the squark hierarchy and the nature of the LSP. 

Similar studies are performed for benchmark Model C. The result for Model C is shown in 
Fig. |4j Similar accuracies of the measurements are obtained with same integrated luminosity. In 
particular, we observe that we should be able to distinguish Model B and C from these results alone 
with about 5 fb _1 of integrated luminosity. 
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As a final note, we add that the even for the correct mass template the ratio Br(g — > tt)/ Br(g — ► 
bb) is slightly underestimated in all cases. This is because the data actually contains a small 
percentage of events in which the gluino decayed to first and second generation quarks. As those 
decays looks most similar the decay into two b-jets, the fit tends to slightly overestimate Br(g — > bb), 
and therefore underestimating Br{g — > tt)/Br(g — ► bb). 
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Figure 5: Results for fits of ratio Br(g — > tb+ib) / Br (g — > tt) for benchmark Model-D. The left panel 
shows the result with correct mass hypothesis. In the right panel, the efficiencies are calculated for 
five different mass templates show in the legend. The solid horizontal line gives the actual values 
of the branching ratios. Errors are la and include errors for subtracting off the tt background. 



The study of gluino decay for benchmark model D is presented in Fig. [5l In this case, we are 
interested in ratio Br(g — > tb + ib) /Br(g — ► tt). For a wino-LSP model like Model D, we expect this 
ratio to be large, which can indeed be experimentally verified with moderate luminosity, as shown 
in the figure. 

We finally consider how well we can chose our mass hypotheses based on available experimental 
data. This is certainly an important issue since significantly wrong mass hypotheses lead to mis- 
leading results. To begin with, we study the dependence of our result on the mass hypothesis in 
more detail. 

As the mass gap between the gluino and LSP is tightened, the events will have a harder time 
satisfying the missing energy cut we imposed. This more significantly affects events of the form 
g — ► tb + fir than g — ► bb + flr, since the tops will use more of the gluinos' energy than the bottoms. 
Thus a tighter mass gap used in our template underestimates the ratio of efficiencies £titt/ e ttbb> see 
Figure [6l The fitter then adjusts for this low efficiency by fitting more tttt events relative ttbb to 
the signature counts 

Br(gg^tttt) _ Br(g -> tt) 2 _ e tM 

Br{gg -» ttbb) ~ 2Br(g -> tt)Br{g -» bb) ~ e tM ' 

Thus, in Models B and C as we increase(decrease) the gluino mass in our templates we tend to 
underestimate(overestimate) the branching ratio Br(g — > tt)/Br{g — > bb). 

We see that the change in the mass gap between gluino and the LSP can account for most of 
the variation in the fit result. For example, extrapolating from our study, if our assumption of 
gluino mass is off by more than 200 GeV, the result for benchmark Model B will indeed look similar 
to that of Model C. However, since we expect to have significant excess in multiple channels, we 
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Figure 6: Dependence of efficiencies on the mass gap Mg — Mlsp to pass the missing energy cut, 
MET > 100 GeV. This dcpcndance is the dominant effect for the variation of fits on the mass 
hypothesis. 

observe that we should also have a significant amount of information about the mass scales, and 
in particular the mass gap M§ — Mlsp of the new physics particles already at early stages of 
LHC running. For example, we would expect simple transverse variables in channels with multiple 
leptons (and hence lower Standard Model background) should already provide some indication of 
the mass scales involved. 

Indeed, visible differences can be seen between different gluino masses when we histogram the 
effective mass, M e s, defined as the scalar sum of the transverse momentum for all objects in the 
an event 

M eff = J2p T ,r+^T- 

i 

To demonstrate this, we plotted the effective mass for gg — > tttt for the 5 gluino masses used in 
fits above, in Figure [7] The histogram curves move to lower energies (right to left) as the gluino 
mass is lowered. There is also a noticeable change in the effective mass spread, which we quantify 
as the location of middle 20% quantile. The bars at the bottom show the location of the inner 20% 
quantile; the highest one is for the largest gluino mass, the second highest for second largest gluino 
mass, etc. The fact that the bars move right to left as the gluino mass is lowered indicates that the 
median effective mass is decreasing, while the fact that the bars are shrinking indicate the effective 
mass has less variation as the gluino mass decreases. 

This analysis was only carried out for gg — > titi, so we still need to demonstrate that the effective 
mass is also independent of the gluino decay. To do this we plotted the effective mass for the four 
models considered in this paper, in Figure [5] The bars at the bottom show the spread, from top to 
bottom, for Model-A, Model-B, Model-C, and Model-D. As can be seen, the four histogram curves 
are remarkably similar, and have no noticeable differences in the median or spread. 

Notice that since our fit yield a set of values for ^JoggC -Br c haimci with any given mass hypothesis. 
It can by itself provide a consistency check on the hypothesis. In particular, we can get a lower limit 
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Figure 7: Effective mass distribution, of gg — > frfi events for the 5 gluino masses used in fits. The 
bars at the bottom show the location of the inner 20% quantile; the highest one is for the largest 
gluino mass, the second highest for second largest gluino mass, etc. The histograms are normalized 
so that the total area is unity. 




Figure 8: Effective mass distribution, of Models-(A,B,C,D). The bars at the bottom show the 
location of the inner 20% quantile. The histograms are normalized so that the total area is unity. 



on the gluino pair production cross section, by summing and squaring the fit values y/aggC -Br c hannei 
Ogg > £ _1 X! (V (TggCBru + y/ a gg CBr bb + y/ 'oggCBr tb ) 2 . 

We can obtain cross sections for gluino pair production in the case of decoupled scalars, and rule 
out some incorrect gluino masses used in the templates. For example, in Model C (see Figure [9|), 
we can rule out the nig = 755 GeV mass template at 200 pb _1 of data at la certainty, and begin 
to rule out the nig = 700 GeV mass template at approximately 1 fb _1 of data at la certainty. If 
the actual gluino is much lower than 600 GeV then there must be a significant branching fraction 
of the gluino that is not contributing the channels we used in our fits, such as gluino decays into 
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first and second generation quarks. 
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Figure 9: Gluino production cross section obtained by summing the individual branching ratios. 
We also show the theoretical cross sections at NNLO for gluino pair production for masses used 
in the templates. Notice that at high integrated luminosity we can begin to rule out heavier mass 
hypothesis 



6 Conclusion 

We have studied the LHC signals of pair produced light gluinos which decay dominantly through 
g — > tt + Xi 9 - * bb + X, an d g — > tb + x- We conclude that an early discovery of new physics in 
this scenario is possible due to significant excesses expected in multi-lepton multi-bottom channels. 
Measuring relative branching ratios of gluino decay into tt, bb and tb channels is essential to extract 
information about the underlying model. The crucial step in such a measurement is identify top 
multiplicity in the signal events. We show that direct reconstruction, while useful in gathering 
evidence for the existence top quark in decay products, is not sufficient to measure the number 
of top quarks in the event effectively. We proposed and studied a method based on fitting a set 
of branching ratios to a collection of experimental observables, most of them inclusive counts. 
Efficiencies for identifying a particular final state resulting from certain underlying decay topology 
are estimated by simulating corresponding templates. We conclude that this method will allow us 
to learn about gluino decay branching ratios with roughly 10 fb -1 of integrated luminosity. We 
verified that, in combination with earlier information on the mass spectrum of gluino and the LSP, 
we can obtain a reliable measurement by choosing appropriate mass hypotheses. We emphasize 
that the main advantage of our method is that it allows us to use a large number of channels, many 
of them with multiple leptons and bottoms, in which we expect to see excesses during the early 
stage of LHC in models of this sort. 

We showed that we can expect to gain enough information about the mass spectrum, in par- 
ticular the mass gap Mg — Mlsp, from simple observables like effective mass and demonstrated 
that, with our fitting method, it is also possible to obtain an estimate of the gluino production 
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cross section which in turn gives us very valuable information of the gluino mass, which also gives 
a consistency check of the assumption we made in our measurement of the branching ratio. 

In the paper, we have considered only benchmark models with only a set of simple decay chains. 
More complicated models will certainly contain channels which requires further study [26j [27] ■ For 
example, a decay chain which contains ...£ — > bC followed by C — > WN has the same set of final 
state particles as top decay. On the other hand, we have also only used counting signatures in 
our fit. Inclusion of more kinematical variables may improve our ability of discerning other decay 
topologies. For example, in the decay of stop mentioned above, the kinematics of b and W will be 
in general different from the case of top decay. 

One key factor entering our lepton efficiency is the isolation requirement. This is particularly 
significant in our case since we expect to have a lot of hadronic activity in these top rich events. We 
have only implemented a simple and commonly used isolation criterion based on hadronic activity 
within a narrow cone around the lepton. However, one can in principle improve and optimize the 
isolation cuts by including additional alternative isolation criteria, such as by requiring the invariant 
mass of lepton and the hadronic activity ma > m C ut- This cut is effective since a lepton from heavy 
flavor decay is typically soft and therefore gives a small m« i; while it is expected to be larger for 
accidental overlap between lepton and jet. 
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